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Abstract 

We propose a general interpretation for long-range correlation effects in the activity 
and volatility of financial markets. This interpretation is based on the fact that the 
choice between 'active' and 'inactive' strategies is subordinated to random-walk like 
processes. We numerically demonstrate our scenario in the framework of simplified 
market models, such as the Minority Game model with an inactive strategy, or a 
more sophisticated version that includes some price dynamics. We show that real 
market data can be surprisingly well accounted for by these simple models. 



1 Introduction 



When looking at the time series of the daily returns in liquids markets, it ap- 
pears to the naked eye that the observed behaviour is markedly non Gaussian: 
fluctuations are intermittent and strong, localized outbursts of volatility can 
be identified. This fact, known as volatility clustering [1-3], can be analyzed 
more quantitatively by looking at the daily volatility o t (denned for example, 
as the average squared high frequency return). This quantity follows a log- 
normal distribution [5], and its temporal correlation function (er t er t+r ) can be 
fitted by an inverse power of the lag r, with a rather small exponent in the 
range 0.1 — 0.3 [1,4-6]. This suggests that there is no characteristic time scale 
for volatility fluctuations: outbursts of market activity can persist for rather 
short times (say a few days), but also for much longer times, months or even 
years. Besides, these long ranged volatility correlations are observed on many 
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different financial markets and over different periods of time, with qualita- 
tively similar features: stocks, currencies, commodities or interest rates (see 
for example the bund contract discussed in [3]). This suggests that a common 
mechanism is at the origin of this rather universal phenomenon. This univer- 
sality, which is also observed for other 'stylized facts' of financial markets (for 
example the power law behaviour of the return distribution) has attracted 
much attention, especially among physicists. It indeed supports the idea of 
the market as a complex interacting system of the kind usually studied in 
physics, where complex behaviour arises from individual actions, not crucially 
depending on the details of the microscopic interactions. 

In this line of thought, the approach of agent based modelling has been in- 
tensively pursued, and important insights into market dynamics have recently 
been gained [8,10-18]. These models postulate some simple behaviour at the 
level of the agents and investigate the resulting price dynamics. Their aim is to 
reproduce real market phenomenology and understand its distinctive features 
in terms of microscopic mechanisms. Of course the possibilities of modelling 
are a priori numerous, since very little is known about the relevant ingredients 
needed to construct a realistic artificial market. In this context two possible 
strategies of investigation can be followed. The first is to look for very simple 
models, which may be in some respects unrealistic, but where both analyti- 
cal resolution and intuitive understanding can be reached. Alternatively, one 
can focus on more complex models, trying to select the necessary microscopic 
structure needed to reproduce the stylized facts, and identify the 'universality 
class' of real markets. 

In this paper we consider models of the two types and we look for the micro- 
scopic origin of volatility clustering. We propose a simple and robust mecha- 
nism to account for the appearance of long-ranged correlations in the above 
simplified models [19]. We then argue that this mechanism also very naturally 
operates in real financial markets, and accounts well for the empirical find- 
ings. Similar ideas (although quantitatively different) were recently discussed 
in [20,15]. 



2 Simple models: the Minority Game 

As an example of the first family of models we consider the Minority Game 
(MG) and its variants [10,21,22,25-27]. This model describes the behaviour 
of competing agents that can choose between different individual strategies 
as a function of their past performance. It was originally introduced to model 
bounded rationality and inductive reasoning behaviour [9,10], but it has since 
then become a paradigm for systems with heterogeneous adaptive agents. 
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The setup of the model is very simple. At a given time each agent % can 
take two possible actions Oj = ±1. Once that the global action A(t) = J2i a i 
is determined, those agents who are in the minority group win whereas the 
others lose. The game is repeated at each time step. It would of course be 
trivial if the agents took their actions by tossing a coin. In fact it is not, since 
the opinion formation mechanism is not trivial at all. Its distinctive properties 
are: 

• Heterogeneities: Each agent i is endowed with s = 1,...,S "strategies" of 
action which are randomly chosen among all possible strategies at the initial 
time, and are kept fixed throughout the game. 

• Information: A public and common information J, (for example the m-steps 
past history [10], or an exogenous signal [25]) is available at each time step. 
Each agent processes this information (the same for all) using one of his 
strategies, and decides of his action. 

• Adaptivity: Agents evaluate the performance of their strategies. Agent i 
gives scores to each of his strategies according to its observed predictive 
power: 

P hS (t+l)=P hS (t)-a hS A(t) , (1) 

where s is a strategy label, and a simple minority payoff function — aA has 
been implemented. Then, at each time t the best strategy, the one with the 
highest score, is actually used to decide the action. 

The collective behaviour of the MG is extremely rich and has been extensively 
studied in numerical and analytical works. The main feature is the possibil- 
ity, by tuning the parameter a given by the ratio of the number of different 
'patterns' of information to the number of agents, to obtain an 'efficient' or 
an 'inefficient' behaviour (in a sense very close to the analogous concept in 
economics). However, in its original version, this model is rather remote from 
financial markets; in particular, there is no price dynamics. Several attempts 
have been made to generalize it and construct more realistic market models 
[21,22]. As first noticed in [22], if one also allows the agents to be inactive, 
intermittent volatility fluctuations can be generated. This fact confirms in the 
context of generalized MG what has been observed also for real markets, that 
is that volatility fluctuations are related to activity correlations [31,30]. For 
what concerns our main point, i.e. to explain volatility clustering, this means 
that we can focus on the temporal pattern of activity in these models. From 
this point of view, the simplest model which exhibits long ranged activity cor- 
relations is the MG with inactive strategies. Here we just consider the original 
setup of the MG, but we give to each agent one more 'inactive' strategy whose 
prediction is always ctj = (see also [23,22,24]). The number of active agents 
will then fluctuate from time to time, and for this reason we shall refer to 
this version of the model as Grand- Canonical MG. As we shall see in the next 
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sections, in this case it is rather intuitive to understand how the microscopic 
individual patterns determine a specific behaviour of the activity correlations. 



3 More complicated market models 

In order to go beyond the above mentioned limits of the MG we have studied a 
more complicated market model that, while retaining the non trivial opinion 
formation structure of the MG, allows traders to switch between a 'bond' 
market (representing the inactive state) and a 'stock' market, and accounts 
properly both for their wealth balance and for market clearing. 
As in the MG each agent has S random active strategies plus an inactive one. 
Besides he owns a number 0j of stocks and Bi of bonds which he can trade 
at each time step. The price X(t) of the stock evolves in time according to a 
supply /demand dynamics which is defined through the following steps: 

• Information: As in the MG at each time t an information is given to the 
agents. We choose it to be given by the m last steps of the past history of the 
return time series: I(t) = (x(t — m), • • • , x(t — 1)) with x(t) = sign(X(t) — 
X(t — 1)). In this sense, our traders are chartists on short time scales, and 
take their decision based on the past pattern of price changes (note that 
this decision is not the same for different agents). 

• Opinion formation mechanism: As in the MG agents use adaptive strategies 
to take a decision, which in this case can be 'buy', 'sell', or 'do nothing'. 
Besides, each agent wishes to buy/sell a quantity proportional to his current 
belongings, that is ai = gBi(t)/X(t) if he wants to buys, Oj = —g4>i(t) if he 
wants to sell. Also, each agent can act as a fundamentalist with a certain 
probability Pf on d that depends on the difference between the observed price 
and some 'reference' price that grows with the risk-free interest rate. This 
leads, on the long run, to a mean reverting behaviour of the price towards 
its reference value. 

• Price formation mechanism: Once the individual offers are made the 
price update is defined by [11,12]: 

5X(t) _ X(t + 1)-X(t) _ _ 

~x(t)- x{t) " A 4^ (t) ' (2) 

where A is a measure of the liquidity of the market. Highly liquid markets 
correspond to small A. 

• Market clearing mechanism: We consider a two steps dynamics. First the de- 
cisions/orders are made and the price is updated. Then a matching between 
supply/offer is realized, which may leave a certain amount of unfulfilled or- 
ders. More precisely, the amount of stocks 5(f)* actually traded by agent i is 
g e Bi(t)/X(t+ 1) if he buys and —g e 4>i{t) if ne sens ; where g e is such that 
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J^i^i = 0- Unfulfilled orders are then removed from the order book. 
Wealth dynamics: After this 'market clearing' the individual wealths are 
updated using the quantity 5(f)* of stocks actually traded by agent i, i.e.: 

<f>i(t + 1) = <j>i(t) + 8$ B t {t + 1) = 5,(0(1 + p) - 5<t>*X{t + 1), (3) 

where p represents the interest rate. The total wealth of agent % is therefore 
w i (t) = <p i (t)X(t) + B l (t). 

Scores update: The scores of the strategies are finally updated according to 
their real profitability over the interest rate p: 



P i ,s(t + l)=iJiP ita (t)-a itS (t-l) 



5X(t) 



(4) 



When p, < 1, only the recent past is used to assess the performance of the 
strategies. 

The phenomenology of this model is very rich and a detailed account of our 
results will be published separately [28]. From the previous description we 
see that there are at least four parameters which enter crucially in the model: 
p,g,X and one entering the precise definition of pf on d- Tuning these parameters 
we observe two qualitatively different regimes: 

1) An Oscillatory Regime (for small values of A and g), where speculative 
bubbles are formed, and finally collapse in sudden crashes induced by the 
fundamentalist behaviour. In this regime, markets are not efficient, and a 
large fraction of the orders is (on average) unfulfilled. 

2) A Turbulent regime (for large A and g) where the 'stylized' facts of liquid 
markets are well reproduced: the market is approximately efficient (although 
some persistent or anti-persistent correlations survive), the returns follow a 
power law distribution, and volatility clustering is present. 

Both these two regimes present interesting features that can be analyzed. 
However, for the objective of the present paper, we focus here on the second 
one. In this regime both the volatility and the activity show power-law like 
correlations (in time) which are quantitatively very similar to those observed 
in the much simpler Grand- Canonical MG: see Figures 1 and 2 below. This 
strongly indicates that volatility/activity correlations are uniquely determined 
by the opinion formation structure which is common to the two models. In 
the next section we propose an explanation of how this can happen. 



4 A simple universal mechanism 



In the above models, scores are attributed by agents to their possible strate- 
gies, as a function of their past performance. In a region of the parameter 
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space where these models lead to an efficient market, the autocorrelation of 
the price increments is close to zero, which means that to a first approxima- 
tion, no strategy can on average be profitable. This implies that the difference 
between the score of two strategies will locally behave, as a function of time, 
as a random walk. Now, the switch between two strategies occurs when their 
scores cross. Therefore, in the case where each agent has two strategies, say 
one 'active' (trading in the market) and one 'inactive' (holding bonds), the 
survival time of any one of these strategies will be given by the return time of a 
random walk (the difference between the scores of the two strategies) to zero. 
The interesting point is that these return times are well known to be power-law 
distributed (see below): this leads to the non trivial behaviour of the volume 
autocorrelation function. In other words, the very fact that agents compare 
the performance of two strategies on a random signal leads to a multi-time 
scale situation. More formally, let us define the quantity 9i(t) that is equal to 
1 if agent i is active at time t, and if inactive. The total activity is given by 
a(t) = ^2i9i(t). The time autocorrelation of the activity is given by: 

C a (t,t>) = (a(t)a(ti)) = Y,mm'))- (5) 

We will actually use in the following the so-called activity variogram, directly 
related to the autocorrelation through: 

V a (t,t') = ([a(t) - a(t')] 2 ) = C a (t,t) + C a (t',t') - 2C a (M')- (6) 



One can consider two extreme cases which lead to the same result, up to a 
multiplicative constant: (a) agents follow completely different strategies and 
have independent activity patterns, i.e. {9i9j) oc 5ij or (b) agents follow very 
similar strategies, for example by all comparing the performance of stocks to 
that of bonds, in which case Q i = 9j. In both cases, C a (t,t') is proportional 
to (9i(t)9i(t')) . This quantity can be computed in terms of the distribution 
P(s) of the survival time s of the strategies. More details can be found in 
[19,29]. When P(s) has a finite first moment (s) (finite average lifetime of the 
strategies), then C a (t,t') is stationary, i.e. it only depends on the difference 
t — if — t. Introducing the Laplace transforms CC a (E) and CP(E) of C a {r) 
and P(s), the general relation between the two quantities reads [29]: 

For an unconfined random walk, the return time distribution decays as s~ 3 / 2 
for large s and therefore its first moment is infinite. However, in all the mod- 
els mentioned above, there exist 'restoring' forces which effectively confine the 
scores to a finite interval [28]. This can be attributed, both in the case of the 
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Fig. 1. 1-a: Cumulative distributions for the survival times of the active strategy in 
the Grand-Canonical MG with an inactive strategy and two active strategies, for 
five different agents. The parameter a is set to a = 0.51 (a c ~ 1.), and the number 
of agents is 501. 1-b: Volume of activity (number of active agents) as a function 
of time for the MG with the same values of parameters. Inset: The corresponding 
activity variogram as a function of the lag r, in a log-log plot to emphasize the yfr 
singularity at small r's. 

MG or of more realistic market models, to 'market impact', which means that 
good strategies tend to deteriorate because of their very use, or to the finite 
memory time of market operators used to assess their strategies (cf. the pa- 
rameter \i defined above). The consequence of these effects is to truncate the 
s~ 3//2 tail for values of s larger than a certain equilibrium time so- Therefore, 
the first moment of P(s) actually exists, such that Eq. (7) is valid. Neverthe- 
less, one can see from Eq. (7) that the characteristic s~ 3 / 2 behaviour of P(s) 
for short time scales leads to ECC a (E) ~ 1 - B/\/E + ... for s 1 < E < 1. 
This in turn leads to a singular behaviour for the variogram V a (r) at small 
r's, as V a (r) oc y/r, before saturating to a finite value for r ~ s . Intuitively, 
this means that the probability for the activity to have changed significantly 
between t and t + r is proportional to J T ds sP(s) oc yfr (for r <^ s ), where 
sP(s) is the probability to be at time t playing a strategy with lifetime s. 

Let us illustrate and justify this general scenario for the MG with an inactive 
strategy. In this model [23,28] there is a critical value a c of the relevant param- 
eter above which all agents finally become inactive. However, below this value, 
the system is always 'efficient' and our assumption about the local random- 
walk behaviour of the strategy scores is thus well grounded. However, in order 
to confirm this point we show in Figure 1-a the cumulative distribution of 
the survival time of the active strategy, for five different agents. In the case 
of a pure random-walk this distribution should behave as s^ 1 / 2 . We can see 
that this is precisely what happens up to a certain time scale s (which does 
not fluctuate much from agent to agent) above which the distributions are 
truncated. In Figure 1-b, for a value of a smaller than a c , we have plotted 
the volume of activity as a function of the time. In the inset, the activity var- 
iogram V a (r) reveals the characteristic y/r singularity discussed above, before 
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saturating for large r (~ sq). This sfr singularity is present in the whole active 
phase a < a c ; s is found to be proportional to N (for ji = 1). The analogous 
variograms for the volume of activity and for the volatility of the more realistic 
market model of section 3 are shown in Figure 2, together with the result of 
the MG. The two behaviors are very similar and confirm the universality of 
this result [28]. Note that the activity variogram of the MG model can be very 
accurately fitted by the following simple form: 



We note that a similar mechanism might also be present in other models 
where agents can switch between different strategies or classes of strategies. 
For example, in [15] each agent can behave either as a fundamentalist or as 
a trend follower, switching between the two strateg 1GS clS db function of their 
relative performance. For this model it has been observed that the activity 
bursts are indeed associated to a large number of agents switching from one 
behaviour to the other. The importance of the fact that these strategies are 
on average equivalent was also clearly stressed [15] (see also [20]). It would be 
interesting to see if the activity variogram in this model can also be fitted using 
Eq. (8). Another example is the adaptive evolutionary system of [16], where 
agents can use different expectation functions to estimate future prices and 
dividends (and then act differently according to a mean- variance optimization 
procedure) and choose between them according to a evolutionary performance 
measure. Q 

It is interesting to compare the above results with real market data. Figure 
3 shows the volume of activity on the S&P 500 futures contract in the years 
1985-1998. On the same graph, we have reproduced the MG fit Eq. (8). Both 
the time scale and the volume scale (arbitrary in the MG model) have been 
adjusted to get the best agreement. Furthermore, a constant has been added 
to V a (corresponding to a <5 Tj o contribution to C a ), to account for the fact that 
part of the trading activity is certainly white noise (e.g. motivated by news, 
or by other non strategic causes). As can be seen, the agreement is rather 
good. Most significant is the clear yfr behaviour at small r (r < 50 days). We 
therefore suggest that the effect captured by the MG or more sophisticated 
variants (Figures 1 and 2), namely the subordination of the activity on random 
walk like signals, is also present in real markets. It seems to us that this makes 
perfect sense since market participants indeed compare the results of different 

1 In this model the presence of positive dividends on stocks is crucial to implement 
a reasonable market dynamics. In the market model described above, on the other 
hand, we do not expect dividends to change substantially the dynamics since a 
dividend at time t is in fact followed by an immediate decrease of the price X(t) of 
the same amount, leaving the scores updates unchanged. 
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Fig. 2. Returns (relative price differences) as a function of time in the market model 
described in the text (section 3). Large frame: The corresponding volatility vari- 
ogram as a function of the lag r, in a log- log plot to emphasize the \fr singularity 
at small r's. The full line is the fit Eq. (8) 
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Fig. 3. Total daily volume of activity (number of trades) on the S&P 500 futures 
contracts in the years 1985-1998: compare to Figure 1 and 2. Inset: Corresponding 
variogram (diamonds) as a function of the square-root of the lag. Note the clear 
linear behaviour for small y/r. The full line is the MG fit Eq. (8), with both axis 
rescaled and a constant added to account for the presence of 'white noise' trading. 
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Fig. 4. Variogram of the log-volatility, (log 2 (cj t /cr t+T )) as a function of r, aver- 
aged over 17 different stock indices (American, European, Asian). The upper panel 
shows a fit with a power law of time, the middle panel a fit with the logarithmic 
(multi- fractal) prediction, whereas the lower panel shows the MG fit Eq. (8), with 
again both axis rescaled and a constant added to account for the presence of 'white 
noise' trading. Note that the three models lead to nearly indistinguishable fits. 

strategies to decide whether they should remain active in a market or leave 
it. Note furthermore that although our scenario is based on the comparison 
between the scores of strategies, similar results would be obtained if the volume 
was subordinated to the difference between the price and certain 'psychological 
levels' (i.e. the value 1000 for the S&P, etc.). It is indeed reasonable that such 
levels play a role in determining the activity on financial markets (see however 
[32] for a discussion of this point). 

Since the volatility and the volume of activity are strongly correlated also 
in financial markets [30,31], our interpretation should carry over to volatility 
fluctuations as well. This is illustrated in Figure 4 (lower panel), where the 
variogram of the log- volatility for major stock indices is shown, together with 
the very same MG fit. Again, the agreement is very good. We have however 
also shown for comparison the prediction of the multi-fractal model of ref. [6] 
V a (r) = 2A 2 log(r/ro) (middle panel). We note that the two models, although 
very different, lead to nearly indistinguishable numerical fits over the time 
scale considered. In the literature, it is also customary to fit the volatility cor- 
relations with power laws. We have therefore also shown in the upper panel the 
corresponding fit for the variogram V a (r) = c — Ci/t u , with v = 0.1. This fit 
is also consistent with the data (and is numerically very close to a logarithmic 
function because v is small). However, neither the power-law nor the logarithm 
have clear microscopic motivations. If, as we believe, the simple mechanism 
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that we have advocated is at the origin of the correlations behaviour then the 
truly universal feature is in the \fr behaviour of the variogram for small r, 
and not in the long time behaviour of the correlation function. In other terms, 
the exponent v might be an effective, non universal exponent masking the true 
universal 1/2 exponent that describes the initial increase of the variogram. 

In summary, we have proposed a very general interpretation for long-range 
correlation effects in the activity and volatility of financial markets. This in- 
terpretation is based on the fact that the choice between different strategies 
is subordinated to random-walk like processes. We have numerically demon- 
strated our scenario in the framework of simplified market models, and showed 
that, somewhat surprisingly, real market data can actually be quite accurately 
accounted for by these simple models (see Figs 3 and 4). 
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